General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



E85'1003y 




CONTINENTAL AND OCEANIC CRUSTAL 
MAGNETIZATION MODELLING 


A SEMIANNUAL STATUS REPORT 


(E85-1UU34 «AbA-CE- 1741^8) CC fall iiENXAL AhC 
CCi;;A.NIC CEU8IAL HAGfaElIZATlCfa ECDELLING 
SemidiiMiuai 8tdtu;:i Qepcct (axaiii Dniv.) 11 p 
ilC A02/ilr A01 CStL 08C 

Gi/** 3 


C. G. A. Harrison and K. L. Hayling 


Roaenstlel School of Marine and Ataospherlc Science 
University of Miami 
4600 Rlckenbacker Causeway 
Miami, Florida 33149 


fad5-i;413 

UncidS 

00034 


1 


Inversion of magnetic data from the MAGSAT satellite, to arrive at 
Intensities of magnetization of the Earth's crust, has been performed by 
two different methods. 

The first method uses a spherical harmonic model of the magnetic 
field. The coefficients believed to represent sources In the Earth's 
crust can then be Inverted to arrive at vertical dipole moments per unit 
area at the Earth's surface (Chapman and Bartels, 1940). The spherical 
harmonic models we have used contain coefficients of degrees of harmonics 
up to 23 (Langel and Estes, 1982) and 29 (Cain et al., 1984). The dipole 
moment per unit area for a surface element can then be determined by 
suBBilng the contribution for each Individual degree of harmonic. The 
magnetic moments have been calculated for continental and oceanic areas 
separately as well as over certain latitudinal segments. 

We wanted to Investigate If there are any differences between 
continental and oceanic areas, which can be suspected from the great 
differences between these In terms of crustal thickness and composition. 
The difference between high and low latitudes should differ by 
approximately a factor of 2, because the strength of the geomagnetic field 
increases by this factor (30,000 to 65,000 nT) when going from equator to 
pole. Table 1 shows the results for the Langel and Estes and the Cain et 
al. models. Neither one of these shows any greater differences between 
continental and oceanic magnetization, and no latitudinal variation can be 
observed. In the case of the relation between oceanic and continental 
magnetization, the unit ratio would mean that the greater thickness of Che 
continental crust, 30-40 km, compared with the oceanic crust, 6 km, is 
compensated by the same but inverse factor of susceptibility of 


continental compared with oceanic rocks* There la no reason why the 
magnetic field, believed to be generated In the core or at the core-mantle 
boundary, should have any preference for continental or oceanic crust In 
terms of Intensity. Even more surprising Is the lack of variation for 
different latitudes* This means that the crust Is thinner and/or has a 
lower average susceptibility towards the poles* 

The core magnetic field Is believed to be represented by degrees of 
harmonics 1 through 13, of a spherical harmonic model* A pronounced 
change In slope In the power spectrun of such models, using the expression 
from Mauersberger (1956) and Lowes (1966), Is occurring between degrees 13 
and lA (Langel and Estes, 1982)* The other main reason for this choice Is 
the white spectrum this band width represents at the core mantle boundary 
In terms of energy (Langel and Estes, 1982)* Figure 1 shows that the 
energy, calculated by the Lowes function at the core mantle-boundary, has 
been constant during the last 2000 years If degrees of harmonics 1 through 
18 are Included, Indicating that the core field Is better represented by a 
wider span of harmonic coefficients* 

The next analysis with magnetization Intensities of the Earth's crust 
from spherical harmonic models was made by using narrower ranges of 
degrees of harmonics, thus assuming that degrees higher than 13 are 
present In the core field signal* Figure 2 Illustrates the vertical 
dipole moment per unit aea for oceanic and continental areas when 
different spherical harmonic windows, as well as the latitudinal 
difference, have been Included* A closer correlation with the theoretical 
variation as a function of latitude la obtained when spherical harmonic 
models with coefficients 17 or higher are representing the crustal 
magnetization* 
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These results sre of the greatest Importance when a data set Is to be 
derived to represent the anomalous field arising from structural and 
compositional variations in the Earth's crust. Such an anomalous magnetic 
field has been used in our second approach to model crustal 
magnetizations, in which inversion is performed by applying the equivalent 
aource technique (Dampney, 1969). The inversions have been performed by 
slightly modified Gauss- Jordan eliminations. We have uaed the scalar 
anomalous magnetic field of Langel et al. (1982) to investigate an entire 
ocean basin, the North and Equatorial Atlantic. The magnetization models 
have been adjusted with a certain amount of annihllator (Parker and 
Heustis, 1974) to arrive at entirely positive or zero magnetizations, a 
necessity if induced magnetization is the leading source for crustal 
magnetization. Remanent magnetization would also Justify adjustment with 
an annihllator, when the magnetic reversal history of the area is 
considered. Local areas of the region have been inverted separately 
because we wanted to investigate local magnetic anomalies and specific 
physiographic features separately, and because the optimal solution (Bott 

m 

and Hutton, 1970) varies from a source spacing of leas than 3 to more than 
4 degrees for different areas. In order to match adjacent magnetization 
models, the annihllator was used so as to give identical or almost 
identical source elements the same intensity of magnetization and entirely 
positive magnetizations. 

The combined magnetization in the North and Equatorial Atlantic is 
shown In Figure 3. Some correlation with geologic structures can be 
recognized, for example the Puerto Rico Trench, the Rio Grande Rise, and 
perhaps the Walvis Ridge. Thus, the lack of anomalies over structures 
like the Bermuda Rise and the Canary and Cape Verde Islands Is surprising. 


4 


During the spreading history of the North Atlantic the Earth's 
magnetic field has been stable during 2 periods of time, with such 
durations that the remanent Imprint should be detectable at satellite 
altitude. Remanent magnetization Is shoim, by direct modelling, to be a 
possible source for long wavelength magnetic anomalies observed over 
oceanic areas when the geometry and location of the quiet zones make them 
observable. The magnetic field from two-dimensional magnetized slabs have 
been calculated to model the quiet zones In the North Atlantic. Figure 4 
shows the theoretical field of a magnetized slab, representing the 
Cretaceous quiet zone (118-83 my, Harland et al., 1983) at 40 degrees 
north. The overall shape of the anomaly can be recognized In the HAGSAT 
scalar anomalous data set, but the Intensity of magnetization of the 
oceanic crust Is on the large side, when considering the results from 
direct studies of oceanic rock samples (e.g. Kent et al., 1978; Dunlop and 
Prevot, 1983). This discrepancy Is even more pronounced In the case of 
Induced magnetization. 

The HAGSAT data, collected at an average altitude of 404 km, should 
enable resolution of signals with wavelengths down to 250 km (Sailor et 
* al., 1982). To Improve resolution, we have started to apply constrained 
least squares Inversion routines to a scalar anomalous data set that only 
permits certain Intensities of mrignetlzatlon In the solution. So far, the 
Simplex algorithm has been used. This method offers the possibility to 
set different constraints for different areas In the same Inversion, but 
has the disadvantage of considerably Increasing the matrices that are to 
be Inverted. Other routines, for example constraining by weights, have 


also been considered 
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Table 1 

Vertical Dipole Kcrient/Area 
OcA) 


Field Representation 



Langel and Eatea, 1982^ 

Cain e^ 

1984^ 

Latitude Band Cont. Ocean Ratio* Cont. 

Ocean 

Ratio* 


• 90-53 

8.63 

9.48 

0.910 

7.36 

7.47 

0.985 

53-37 

5.99 

6.03 

0.993 

6.61 

6.84 

0.966 

37-23 

6.35 

4.90 

1.296 

6.45 

6.60 

0.978 

23-13 

5.96 

6.02 

0.990 

6.45 

6.62 

0.975 

13-0 

5.1b 

6.45 

0.893 

6.37 

6.85 

0.931 

Overall ratio 



1.077^ 



0.990^ 


*Contlnent/Ocean ratio 
'^’Degree 14-23 
^Degree 14-29 
^Areally weighted tneans 


7 


r" 


Figure 


/ 






600 - 

400- 



>- 

C9 

cr 

UJ 

z 

LiJ 


100 - 

90- 

80- 

70- 

60- 

50-. 


40- 


l 



I ‘ ‘ ^ ■ ■ > < ^ 

10 20 


DEGREE OF HARMONIC 


1* Total cumulative energy external to the core plotted aa a 
function of the maximum degree of harmonic uaed In the 
calculation. The lover curve la for the preaent day field ualng 
apherlcal harmonica from Langel and Eatea (1982). The upper cui^e 
la for the field about 2000 yeara ago. during which time the dipole 
component waa about 50Z larger than today. The energy from the 
dipole component la therefore 2.25 tlmea Ita preaent value. The 
harmonica between 2 and 10 were aaaumed to have Increaaed by the 
aame percentage aa the dipole term, becauae the aecular variation, 
aa far aa we know, waa about the aame 2000 yeara ago aa It la 
today. In order to achieve an energy which today la equivalent to 
that In the flrat ten degreea of harmonic 2000 yeara ago. we need 
to uae all harmonica up to and Including 18. 
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Figure 2. Vertical dipole ■oaent per unit area (noraallaed to that In the 
loweat latitudinal band) for (a) continental areas, (b) oceanic 
areaa, and (c) the Earth aa a whole. The latitudinal banda are 
shown along the abscissa. The values expected froa a dipole field 
Intenalty are shown by asterisks. The other syabols show results 
froa different coablnatlons of spherical haraonlcs. -f, degrees 14- 
23; o, degrees 15-23; x, degrees i7-23; □ , degrees 19-23. The 
spherical haraonlc aodel of Langel and Estes (1982) Is used. 
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Plrrure 4. Theoretical Magnetic field froa a aagnetlzed alab at 400 ka 
altitude. The alab la aagnetlzed to 1 A/a, haa a thlcknaaa of 6 
ka, a vldth of 700 ka, and la striking N40* E to the direction of 
the aaln field. 
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